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^SJ ' Abstract 

^^ ' We obtain lower and upper bounds on the heat kernel and Green func- 

vy—j , tions of the Schrodinger operator in a random Gaussian magnetic field 

(^ • and a fixed scalar potential. We apply stochastic Feynman-Kac represen- 

^h ' tation, diamagnetic upper bounds and the Jensen inequality for the lower 

("^ , bound. We show that if the covariance of the electromagnetic (vector) po- 

^^3 ' tential is increasing at large distances then the lower bound is decreasing 

^Z^ ' exponentially fast for large distances and a large time. 
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Introduction 



j^ ' Random magnetic fields appear in many models of physical interest. In the 

^ , fundamental quantum theory the random magnetic field can be considered as a 

^ ' thermal (high temperature) part of the quantum electromagnetic field at finite 

temperature [1][2]. The classical random magnetic field is discussed in optics 
[3] and in a description of an interaction of light with atoms [2]. A random 
magnetic field can arise as a Lagrange multiplier in models of interacting quan- 
ta ' tum particles creating four-fermion (or four-boson)interactions [4] [5]. In the 

Ginzburg-Landau theory of superconductivity [6] when fluctuations are taken 
into account then the electromagnetic field becomes a random Gaussian variable 

[7]. 

The effect of a random electric field in models of quantum mechanics has 
been well elaborated. Anderson localization [8] has been proved for a large 
class of models [9] . It seems that much less is known rigorously concerning the 
localization properties of the Hamiltonians with a random magnetic field (see 
[10] for a recent review; a special case of a magnetic field orthogonal to a plane 
and varying inside the plane is discussed in [11] [12]). Some aspects of localization 
(e.g., the integrated density of states [11]-[12]) can be studied by means of 
the heat kernel of the Schrodinger Hamiltonian in a random electromagnetic 



field. The simplest estimate on the heat kernel follows from the diamagnetic 
inequality [13] [14] which bounds the heat kernel in a magnetic field by that 
without the magnetic field. Such estimates are not interesting if we wish to 
detect the impact of the magnetic field on physical systems. Stronger upper 
bounds on the heat kernel in a (deterministic) magnetic field have been discussed 
in [15] [16] [17] [18] [19]. In these estimates the contribution of the lowest positive 
eigenvalues and eigenfunctions has been estimated. 

We discuss a lower bound on the expectation value of the heat kernel. The 
heat kernel is gauge dependent. We explain which properties of the heat kernel 
do not depend on the choice of the gauge. Our result admits a fast decay of the 
heat kernel for large distances and a large time if the random vector potential has 
growing correlation functions. The well-known example of a constant magnetic 
field [19] shows that the exponential decay is really possible. The effect of 
the magnetic field upon the upper bounds is hard to derive. The classical 
Cwickel-Lieb-Rosenbljum bound on the number of eigenvalues [14] depends on 
momentum and the vector potential A in the combination ]p + A). Hence, 
the dependence on the vector potential A drops out. A more precise bound 
has been derived in [20]-[21] which shows a dependence on the magnetic field. 
Results discussed in [11] - [15] display the discrete spectrum in the asymptotic 
behaviour of the heat kernel for large time and large distances. Although our 
lower bound does not give an exact behaviour of the heat kernel, we nevertheless 
have a feeling that the exponentially decreasing lower bounds reflect an intrinsic 
property of the growing vector potential. We discuss in the last section a simple 
example of a non-trivial bound from above which supports our argument that 
growing vector potentials improve localization. 

The heat kernel defines Green functions in Euclidean quantum field theory 
of an electromagnetic field interacting with a complex scalar field. The relation 
follows from Schwinger's proper time formalism [22]: the correlation functions 
of the scalar field can be calculated from the heat kernel by an integral over the 
(proper) time. We consider the Ginzburg-Landau model with the Lagrangian 



^{ihv + A)^{-iny + A)0 + im^i^p ^ ^|^|4 ^ 1^ 



where Fjk = djAk — dkAj. 

The electromagnetic Lagrangian leads to ultraviolet problems which should 
be irrelevant for a large distance behaviour. We regularize the electromagnetic 
potentials and discuss Euclidean scalar fields in a regular Gaussian electro- 
magnetic field ( the regularity can be achieved by an addition of higher order 
derivatives of the electromagnetic field to the Lagrangian). We obtain lower 
and upper bounds on the correlation functions of the scalar field. Decay of 
correlation functions of the scalar field describes a disappearance of the long 
range correlations of the order parameter in the Ginzburg-Landau model of 
superconductivity [6] [7] . 



2 The imaginary time evolution 

We can apply simple inequalities only for the imaginary time evolution 

|-^ = Ai^ (1) 

where ~1iA{A., V) is the Hamiltonian with a random vector potential A and a 
scalar potential y(the scalar potential could have been random but its eventual 
randomness would not change our results concerning random magnetic fields, 
so we treat it as deterministic) 

nA=-^{-in\7 + Af -V (2) 

2™ 

The imaginary time can be treated as the inverse temperature /3 of the quantum 
Gibbs distribution 

In this section we assume that the scalar potential V is continuous and bounded 
from below and the sample functions of the vector potential A are continuously 
differentiable (almost surely). In the next section we relax the assumptions on 
A (admitting continuous transverse vector potentials). 

Under our assumptions on the potentials for almost every sample path A 
the solution of eq.(l) defines a random semigroup T^ [23] expressed in the form 
[14] [24] 



M^) = (T,^)(x) 



where 



c^P ( I /o^ ^(^ + ^^'s) ° '^^'^s ^ I /o' ^(^ + crbs)ds) V^x + crb^) 



(3) 



(4) 
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and b is the Brownian motion defined as the Gaussian process with the covari- 
ance 

E[hj{s)hk{s')] = 6jkmin{s, s') (5) 

The stochastic integral in eq.(3) is defined in the Stratonovitch sense [29]. We 
have 

a f Aodh{s) = a f Adh{s)+ f ^\/Ads (6) 

where the integral on the rhs is in the Ito sense. Hence, if the vector potential 
is in the transverse gauge then the Stratonovitch and Ito integrals coincide. 

Eq.(3) defines the kernel Kr ' {x, y) as an integrable function on R^ x R^ . 
Under the assumptions that y is a continuous function bounded from below and 
the sample paths of A are continuously differentiable there exists the unique 



kernel continuous on R^ x R^ and determined by the formula [14] [25] (for a 
minimal set of assumptions which allow to derive the stochastic representation 
of the kernel see [26]) 

i^(^^^)(x',x) = cxpM)(x',x) = (27rTa2)-* exp(--i^(x - x')') 

(7) 

E 



exp i- j A(q(s)) o dq^ - - / V{(\s)ds 



here 

q^ =x+(x'-x)-+crVT^a(-) (8) 

T T 

where the Brownian bridge a is the Gaussian process on [0, 1] starting in at 
s = and ending in at s = 1 with the covariancc 

E[aj(s)ak{s')] = 5jks'{l - s) 
if s' < s. It can be expressed by the Brownian motion 

a(.) = (l-.)b(^) 

From the representation (7) there follows the diamagnctic upper bound ( conse- 
quences of the stochastic representation are discussed in [13] [14]; for an analytic 
proof see [28]) 

Proposition 1( Kato inequality) 

Assume that A is a continuously differentiable function and T^ is a continuous 
function bounded from below then the version of the kernel which is continuous 
on R^ X R^ satisfies the inequality 

|i^(^'^)(x,y)|<i^(°^^)(x,y) (9) 

for any t, x, y. 

Let us note that if 

A' = A + Vx (10) 

then 

% i 

exp(T^')(x',x) = exp(-x(x') - -x(x)) exp(r^)(x',x) (11) 

The kernel is not invariant under the gauge transformations. However, 

(V';,exp(T^')V'2) = (V'i,exp(T^)V'2) (12) 

if 

^'j = exp(^x)V', (13) 



We consider a random Gaussian electromagnetic field A' with the mean equal 
zero. We define the covariance G" of A' in an arbitrary gauge as the expectation 
value 

g;,(x,x') = (4(xK(x')) (14) 

It will be convenient to fix the gauge of the random vector potential and sub- 
sequently discuss a dependence of our results on the gauge. From eq.(6) it can 
be seen that the transverse (Landau) gauge divA = is distinguished from the 
point of view of the path integral (if divA. = in a distributional sense then a 
continuity of A is sufficient for the stochastic representation (7) of the kernel). 
If we wish to calculate the heat kernel in another gauge A' then we need the 
gauge transformation with 

X = A^^dwA' (15) 

transforming A' to the transverse gauge. The covariance G of A is related to 
that of A' by the formula 

G,fe(x,x') = (V - 9,a,A-i)(4m - a^a;,A'-i)G;„(x,x') (le) 

We calculate the Gaussian expectation value of the kernel (7) using the definition 
of the Gaussian variable 

(cxp(^(A,f)))^exp(-i((A,f)2)) (17) 

From this formula we can see that a change of the gauge involves the factor 



exp [ -^ I (x (x) - X (x) + / A' o dq 




in the kernel (7). 

3 Bounds on the heat kernel and Green func- 
tions 

We can calculate the expectation value of the kernel over the magnetic field 

(i^(^'^)(x',x)) ^ (cxp(ryl(A,l/))(x',x)) = (2^ra2)-^ exp(--i^(x - x')') 

ZTcr 

£;[exp(-^((^ A(q(,s)) o dq,)2) - i ^ V{cis)ds) 

(19) 
In the formula (19) the Stratonovitch integral can be expressed by the Ito inte- 
gral (the integrals coincide for transverse vector fields). It follows from eq.(19) 
that 

(i^(^^^)(x',x))<i^(°'^)(x',x) 



which is already a consequence of the Kato inequahty (9). 

We rewrite the Ito integral in eq.(19) in the form ( we can prove this equality 
differentiating both sides of eq. (20) below by means of the Ito formula [29] and 
subsequently integrating it again, see also [30]) 



( / Adqf = 2 / Adq, / Adtis' + / A^ — 
Jo Jo Jo Jo ''' 



(20) 



Let us note that 

/ A(q(s)) odcis= / A(q(s))dq, + ^ / divA{q{s))ds 
Jo Jo ^ Jo 

= I r dtvA{ci{s))ds - aV^ r ^A(q(s))b(^) (21) 

^ Jo Jo ^ T — s 

+a^ r A(q(s))(l - -)rfb(^) - r(A(q(s))(x - x')- 

Jo T T-S J„ T 

consists of four terms which could behave in a different way for large distances. 

We shall discuss either a random electromagnetic potential A which is bounded 
(in a certain gauge) or a random electromagnetic potential which is scale invari- 
ant and growing with the distance 

A(Ax) ~ A''A(x) (22) 

where the approximate equality means that the random fields on both sides 
have the same correlation functions. The scale invariance is assumed only for a 
convenience. In general, we can treat the scale invariance (22) as an asymptotic 
behaviour for large distances. 

When 7 > then a scale invariant random field in the transverse gauge (16) 
must have the covariancc (in the Fcynman gauge we could interpret this vector 
field as D-independent Levy's D-dimcnsional Brownian sheets [31]) 

G,fc(x,x') - (<5,,-a,a,A-i)(4r-45;A'-i)(|x|2^ + |xf^-|x-xf^) (23) 

with 7 < 1. The derivatives on the rhs of eq.(23) should be understood in the 
distributional sense. The sample paths of A are Holder continuous with an in- 
dex a arbitrarily close to 7 [27]. A random vector field A which is only Holder 
continuous does not satisfy the conditions for the representation of the kernel 
formulated in sec. 2. However, the kernel (7) for A in the transverse gauge can be 
defined under the weaker condition of Holder continuity [26] . The Ito stochastic 
integral in cq.(7) is defined for any continuous A. The expectation value (7) 
determines a jointly continuous function of r, x, y for continuous potentials A 
and V ( V bounded from below). Then, the kernel (7) defines a random semi- 
group Tt- [23] as a result of the the Markov property of the Brownian motion 



and the multiplicative property of the exponentials [14]. The Hamiltonian —HA 
(together with its domain) is defined as the generator of TV. Alternatively, we 
could introduce a covariance Gn of the vector potential which is infinitely diffcr- 
entiable and depending on a regularizing parameter k . We consider expectation 
values of the kernel with a regularized electromagnetic field. Subsequently, we 
take the limit G^ -^ G oi expectation values (where G is defined in eq.(23)). 
The inequalities discussed in this paper are preserved under such limits. 

We can obtain a lower bound on the heat kernel from the Jensen inequality 
[32]-[33] as applied in cq.(19) to an average over the Brownian motion 

(exp(T^)(x',x)) > (27rra2)-*cxp(-2^(x-x')') 

(24) 
exp ( - ^e[{{J^^ A(q(,s)) o dq.)2)] - i I' V{q{s))ds) 

On the basis of cq.(24) and the diamagnetic inequality (9) (Proposition 1) we 
obtain the following bounds (the terms on the Ihs of eqs.(25) and (28) below 
can be related by some inequalities but we keep this form of the inequalities in 
order to make the origin of these terms visible in such a form as they come from 
eq.(24)) 

Theorem 2 

Assume that c < V < a and the covariance of the electromagnetic field 
(j4j(x)j4fe(x')) in the transverse gauge is continuous and bounded. Then, there 
exists a constant C > and positive constants aj such that 

C(27rTa2)"fexp(-^(x-x')2) 

exp ( - aih^^{x - x.'f - a2h^^x - x'| Vr - aTi^^) (25) 

< (i^r(x',x)) <exp(-|T)(27rr(T2)-#cxp(-^— 2(x-x')^) 



Theorem 3 

i)Assume that A is transverse and continuous and V is continuous and 
bounded from below then the kernel satisfies the upper bound 

(i^(^^^)(x',x)) 

< (27rra2)-* exp(-^(x - x')^) f ^E[exp ( - Jv-(q.))" 

= (27rcr2T)""?exp( (x-x')^) 

2t 

r^ r 2 

ds dy(27r)"*exp(-y)exp(^- ^y(x + (x' - x)s + VTcrs(l - s)y) 

(26) 
ii) Assume that the vector potential in the transverse gauge is scale invariant 
and growing with the scale index 7 (eq.(23)), the scalar potential is a continuous 
function bounded from below and for certain a > and B > 

F(x)<B|x|2^ + a (27) 

then there exists a constant C > and some positive constants aj such that 

(/^.(x',x)) > C(27raV)-f exp(-2^(x-x')2) 

exp ( - aih-\\x\^^ + |xf'')(x - x')2 - a2?i"^ iM^^ + |x'|2-r)|x - x'| V? 

-a3?i-i(|x|2-r + Ix'pTjr - aih-i+''T^+^'\x - x'| - a5h-^+''T^+'^ 

-a6?i-'+^|x - xf T^ - a7?i-i(|x|2/3 + |x'|2/3)t - ash-^+^'r^+f" - ah-\^ 

(28) 
Remarks: 

l.The upper bound (26) does not depend on the vector potential as a conse- 
quence of the Kato inequality (9). The final inequality in eq.(26) follows from 
the Jensen inequality as applied to the time integral. 

2. By a scale transformation (for y = 0) we could obtain (an expectation 
value over a scale invariant magnetic field (23) ) 

(M^'0)(x',x)) ^ (27rcrV)-^exp(--V(x-x')^)cxp(-Ti+'^i^(T-3x,r-3x')) 

(29) 
where the function F has to be determined by an explicit calculation. The lower 
bound (28) is in agreement with the scaling (29); it gives an upper bound for 
the function F. 



The heat kernel will depend on the gauge. However, its diagonal 

(exp(Tyi)(x,x)) >C{2ttt)-^ 

exp ( - 2a3?i~i|x|27r - a5h-^+'>T^+^ - 2a7n-^\x\^^T - asTi-^+^r^+f^ - afi-^r 

(30) 
is gauge invariant. The diagonal of the heat kernel is equal to the Laplace 
transform of the integrated density of states [25] [11]. The integral over the 
diagonal 

(Tr(expM))) = (^ exp (-Te„ (A, V))) = / dx(exp(r^)(x, x)) (31) 

n •' 

is expressing the sum over eigenvalues e„(A, y) of the Haniiltonian —fiA. We 
have 

Corollary 4 

For any r > there exists a constant C > such that 



< (Tr(exp(Tyl))) < (27rTcr2)--T Jrfxexp(-^T/(x)) 



(32) 



here 

^=y(l + -) (33) 

2 p 

with p = max {-J, /3) 
Remarks: 

l.The factor t^'^ in the lower bound on the Ihs of eq.(32) is non-trivial only 
for a small time; for large time the exponential terms decay much faster. If 
y = then the index i' (33) follows already from the scaling (29) ( there is no 
upper bound in eq.(32) if V^ = 0). 

2. The lower bound (32) is a result of the integration over x in eq.(30). The 
upper bound follows from the bound (26) 

J dx{Ki^'^\x,x.)) < (27rTcr2)-# Jrfx£;[exp(- i/J'T/(x + V7cra,)ds 
= (27rTCT2)-f Jdxexp ( - |V(x)) 

(34) 
The upper bound (34) follows from that for J dx J dx{Kr ' (x, x)) (Kato in- 
equality) and has been derived earlier in [13] as a consequence of the Golden- 
Thompson inequality. 

3. The lower bound of Corollary 4 is suggesting that a growing random vector 
field has a similar effect as a growing scalar potential leading to localized states. 
For the harmonic oscillator (with the oscillation frequency w) Trexp{--TH)) = 
(sinh(^))^^. Hence, an increase of the index u in eq.(33) agrees with the exact 
formula (the index v has been discussed also in [34]). However, the exponential 



decrease in the lower bound on the Ihs of eq.(32) does not reflect the exact large 
time behaviour of the trace of the heat kernel of the Hamiltonian with a scalar 
potential. 

4. We can obtain the general formula for {Kr) from the transverse case trans- 
forming a general potential to the transverse one and subsequently calculating 
the average over the gauge function x ^^ i^i eq.(18). The behaviour for large 
distances would not change substantially in Theorem 2 and Theorem 3 if we 
worked with an arbitrary gauge. We have assumed the transverse gauge in or- 
der to avoid difflculties with differentiability of the potentials. Let us explain the 
problem using as an example the square of the first term on the rhs of eq.(21) 



(35) 



\ ds / ds'E[{divA'{qs)divA'{qs'))] 

"'0 

ds ds' E[djdlGjk{qs,qs')] 
Jo Jo 



If the second order derivatives of G are bounded then the term (35) is bounded 
by ct"^ . However, in eq.(23) (without the projection on the transverse part) the 
second order derivative behaves as jxp'''^^ which is singular for 7 < 1. Then, 
the large distance behaviour will be r^jxp'''^^ — T^^'^\T^2x.\'^'f^'^ in agreement 
with the scaling formula (29). 

We consider now the Green functions Qm defined as solutions of the equation 

{-A+^m^)g^m = S (36) 

The Green function can also be defined as the kernel of the inverse operator in 
the Hilbert space of square integrable functions L^(dx) [35]. Then, 



1 f^ 1 

{-A+ -m'^y^ = / dTCxp( m'^T)exp{TA) 



(37) 



By means of an integration over r of the diamagnetic inequality (26) we obtain 
an upper bound for the Green function in a magnetic field in terms of the Green 
function without the magnetic field 

(g(f-^)(x',x))<(?(„°^^)(x',x) 

Under the assumptions of Theorem 3 we obtain the lower bound 

(g(„^^^)(x,x'))>qx-xr^+'exp(-^|x-xf -(|+m)|x-x'|) (38) 

The lower bound for the random magnetic field with the covariance(23) and the 
scalar potential V (27) follows from cq.(28) by an integration over t. A detailed 
estimate of the behaviour of such integrals as a function of |x— x'| is complicated. 
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Without detailed estimates we can obtain an exponential decay in |x — x'| of the 
lower bound for 5m(x, x') as follows from the first r-independent term in the 
exponential on the rhs of eq.(28). It is not clear whether this exponential decay 
comes solely from the unprecise lower bound or if it is an intrinsic property of 
growing vector potentials. 

The Green function is gauge dependent as follows from eq.(ll). We define 
the diagonal 

5Cf.^)(x,x) - g(°'")(x,x) = j^ drexp(-i™V)(i^(A^^)(x,x) - if (°'°)(x, x)) 

which is gauge independent (see [36] for some estimates close to the diagonal ) . 
From eq.(7) 



|g(f ^^)(x,x) - g(„°^°)(x,x)| < j^ drexp(--mV)(27rr)-^ 

/ A(q) odq\+ f dsViq, 
Ja Jo 



(39) 



Hence, if V is continuous and A is Holder continuous with any index 7 > 
(as in eq.(23)) then the diagonal (39) is finite if D < 3. This is so because the 
T- integrand on the rhs of eq.(39) for a small r is bounded by r"""' 2- and for a 
large r it decays exponentially . In fact, we have / j4(q)dq = J(A(q) — A(x))dq 
(because q is a closed path) and for an estimate of the rhs we can apply eqs.(43)- 
(46) below ( see similar estimates in [37]). 

Applying the inequality |(i^)| < (|i^|), the Schwartz inequality and an esti- 
mate of (-B[(/q^ Adq)^]) in cq.(41) below we obtain for 13 < 3 and V > 

|ei°-°)(x,x) - (g(f ^^)(x,x))| < Ci(m)|x|2" + dm) (40) 

where 2a — max{'^, 2/3) and Cj{ra) > 0. 



4 The estimates 

We discuss now estimates leading to the results of Theorem 2 and Theorem 3. 
The upper bound in eqs.(25)-(26) is an elementary consequence of the formula 
(19) ( or the Kato inequality of the Proposition 1). For the lower bound we 
estimate the expectation value on the rhs of eq.(24). An explicit calculation of 
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the average over the electromagnetic field gives 

Jo ~ Jo ~ 
(x-x')i?[G(q(s),q(s'))](x-x') 

Jo Jo T T T - S T- S' 

+ra^ r d(^)(l - -f J2 ElGnHs), q(s))] 

Jo T^-S ^ j 

„2aVT T- r^iJ[b(^)G(q(.),q(.'))](x-x') 
Jo 'T Jo '^ '^-s 

-2o^ r^ r^i5[b(^)G(q(.),q(.'))](x-x') 
Jo T Jo T T-s 

+V7a r^iJ[(x'-x-V^ab(^)) rdb(^)G(q(s)),q(s'))](l--) 

Jq T T - S Jq T - S' T 

(41) 
Note that from the possible nine terms in cq.(41) there remained only six because 

E[ I Usdh{s)] = (42) 

Jo 

if Us depends on b(s') with s' < s (non-anticipating integrals [29]) 

In order to estimate the lower bound in eq.(24) we need some estimates on 
the Ito integrals. We have [29] 

E[{ f tdhsf] = f E[i^]ds (43) 

F{s)^ r{{s'Ms'))dh{s') (44) 

Jo 



Let 



then from the Schwartz inequality 

\E[f dsh{s,h{s))F{s)]\^ < f dsE[h{sf] [ dsE[F{sf] (45) 

Jo Jo Jo 

and for random f and h 

\E[{ r dsh{s, h{s))F{s)])\^ < r ds{E[h{sn{ r dsE[F{sn (46) 

Jo Jo Jo 
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where from eq.(43) 

E[F(sr]= ri(s'Ms')fds' (47) 

Jo 

For the lower bound in eq.(24) we bound each of the six terms in eq.(41) 
from above. Then, we obtain the lower bound for the heat kernel inserting the 
upper bound with the minus sign for each term in the exponential in eq.(24). 

Proof of Theorem 2: 

The upper bound in Theorem 2 is a direct consequence of the Kato inequal- 
ity (Proposition 1). For the lower bound (25) it is sufficient to insert the upper 
bound for each term in eq.(41). The terms without the stochastic integrals in 
eq.(41) can be estimated by means of the Schwartz inequality using the bound- 
edness of V and G whereas for the term with the stochastic integral we apply 
the inequalities (45) and (43). 

Proof of Theorem 3: 

We could obtain the lower bound (28) estimating cq.(41). In particular, 
we could apply the inequality (45) directly to the last term in eq.(41) (then h 
and f in eqs.(43)-(45) do not depend on the electromagnetic field). However, 
it is instructive to return to eq.(21) in order to estimate the product of the 
terms in the square of J Adq directly. Then, h and / depend linearly on the 
electromagnetic field. 

We do not estimate all the six terms which come from the square of J Adq 
in eq.(20)( divA. = 0) but concentrate on three typical terms. The remaining 
three terms can be estimated in a similar way. First, let us consider the last 
term on the rhs of eq.(20) which comes from the square of the stochastic integral 
(20) (it corresponds to the third term on the rhs of cq.(41)) 



ds{E[{A(ci(s))A{ci{s)))] 
= J^ dsJ2E[G,,{q{s)),q{s))] (43) 

Jo 
In eq.(48) the Holder inequality 

la + bpT <22T-i(|a|2'>' + |b|2'^) (49) 

has been applied to the covariance G of eq.(23) with q(s) defined in eq.(8). The 
square of the last term of cq.(21) can be estimated as follows 
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r d. r ^^[^A(q(^))(^, _ ^)A(q(s'))(x' - x))] 
Jo ^ Jo ^ 

Jo 7" Jo ^ 



(50) 



'0 ' Jo 

<ai(|x|2')' + |xf'')(x-x')2 



aelx-x'pcr't'r'' 



where for an estimate of G(q(s)), q(s')) the same set of inequahties has been 
apphed as in eq.(48). 

Next, let us consider a term in the square of J Arfq which is of the form of 
the expression appearing on the Ihs of cq.(46) (cf. cqs.(20)-(21) to see where 
this term comes from) 



/ = crMEilo $A(q(,s))b(^) J^ A(q(.'))(l - t)^^^)]) 



Now 



and 



Ms,b(,s)) = A(q(,s))b( ) 

T — S 



i s\h 



A (q (.'))(! 



Therefore in cq.(46) 



{E[h^]) = E[h{ )G(q(.)),q(.))b(^)] 

T — S T — S 



(51) 
(52) 

(53) 
(54) 



and 



ds / ds' ( E 
Jo 



fMs',b 



Cds rd.'(l- -)2^i?[G,,(q(,s'),q(.'))] 
Jo Jo ^ , 



(55) 



Hence, from eqs.(45)-(46) and the inequality (49) we obtain an estimate on / 

|/|<a3T(|x|2-' + |x'|2^) + a5Ti+^ (56) 

On the basis of eq.(41) and the estimates (48)- (56) it is clear how the lower 
bounds in eqs.(25) and (28) come out ( the estimate in the lower bound (28) on 
the potential V (27) is a simple consequence of the inequalities (49) and (24)). 
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5 Discussion and Outlook 

First of all, let us point out that a deterministic linearly rising vector potential 
really can lead to exponentially decaying Green functions. 

The heat kernel of the Hamiltonian with a constant magnetic field B in 
D = 3 satisfies the inequality [19] [14] (transverse gauge, the z axis is in the 
direction of B and x = (x, y, z)) 

\Kr{x,y,z;x',y',z')\ = {2nT)-^ B{2tt sinhB^)-^ 

I exp ( - ^B^{z - z'f - f coth(Bf )((x - x'f + {y - y'f) + f (V - x'y)) 

< {2TTT)-i B{4:n smh{B^))-^ 

exp ( - ^B^z - z'f -^{{x- x'f + {y - y'f)) 

as 

Bt 
coth( — ) > 1 

By an integration over r (37) we obtain an upper bound on the Green 
function 

\Qm{x,y,z;x',y',z')\ 

< Cexp ( - ^{{x - x'f + {y- y'f) -BI\z- z'\ - m|x - x'|) ^^'^' 

( in the r-integral the inequality sinh -3^ < 5 exp(-^) has been applied). 

The decay of the Green function (57) supports a heuristic argument that 
the term A^ in the Hamiltonian (2) is acting like a potential (see [41] and the 
precise results of [20] [21]). Note that the diagonal 

Kr{x,y,z;x,y,z) ^ F{t) (58) 

being independent of any spatial coordinate is not an integrablc function in any 
of the components of x. 

From the formula for the heat kernel in terms of eigenfunctions and eigen- 
values one can study their dependence on the random magnetic field. For this 
purpose we would need an upper bound for the heat kernel which is stronger 
than the diamagnetic inequality (26) of Theorem 3. We are unable to derive 
such estimates in general. In order to study the localization effects of a random 
magnetic field we investigate a particular model. Let us assume that the mag- 
netic field depends only on coordinates (x, y) of the XY plane. Then, we can 
choose A = (0, 0, ^3(2;, y)). In such a case the Hamiltonian —hA reads 

nA = - — {p + A^f - Vsiz) + —A^y - V2{x, y) (59) 

where A^y is the two-dimensional Laplacian. In eq.(59) we added a potential 
V3(z) ensuring a localization in z. We investigate the conditions on V2 which 
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imply a finite trace of (exp tA) . Wc apply the Golden-Thompson inequality 
[38] [39] ( for a precise formulation and the assumptions see [40]) 



Tr{exp{TAj} < Tr{ exp(^S) exp(TC) cxp(^S) 



(60) 



a A ^ B + C. The rhs of the inequality (32) is a consequence of the Golden- 
Thompson inequality. Now, we choose 



hC 



1 

4to 



and 



We have 



hB- 



ip + A,f-V3{z)-V2ix,y) 



4m 2m 



exp(T,B)(x;x') = (27rcr^T) "'^ exp(- 



1 



/ 



(ipexp(— p(z' — z))E 
n 



exp 



2TCr2 

1 



{x-x'Y' 



1 



2tct2 
{p + A,{q{s))fds 



iv-y'?) 



Ahm Jq 
here q = (gi, 52) is two dimensional (the components defined in eq.(8)) 



(61) 
(62) 



(63) 



exp{TC){x,y,z;x',y',z') = (ttctV) 2 exp{-:^{z - z'f)6{x - x')6{y^- y') 
expi-^V2ix,y))exp{j-{z' - z)A3{x,y))E\cxp ( - j- J^ V3{q3is))ds) 



(64) 



where 



S T S 

93(5) = z + (z' - z)- + \ -aasi-) 
T V 2 r 



From the Golden-Thompson inequality 

{Tr{exp{TA))) < /dx /dx'(cxp(TS)(x,x')exp(TC)(x',x)) (65) 



The expectation value over the magnetic field on the rhs of eq.(65) can explicitly 
be calculated. We perform the calculations in a special case when the potential 
V3 is quadratic 

V3{z)^muj^z^ (66) 

In such a case the expectation value over q^ gives the heat kernel of the harmonic 
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oscillator. After a calculation of integrals over z and z' we obtain 

(rr(cxp(Tyl))) < (27rrcr2)-^(sinli(wr))"i / dxdy exp(-^y2(a;, y)) 

Jdp{exp i^-^k:^ sinh (ujt) (cosh (wr) + 1)"^ (p + A3 (a;, y))^) 

{p + A,{q{.s))fds] 



E 



exp 



Ahm 



< (27rrCT^) 2(sinh(a)T)) ^ dxdyexp{~-V2ix,y)) 



dp l^^{E 



exp 



1 



2mhuj 



sinh(wT)(cosh(u;T) + 1) ^(p + A3,{x, y)f 



Ahm 



{p + A,{ci{s))Y\ ) 



(67) 

The expectation value over the magnetic field on the rhs of eq. (67) can be 
calculated with the result (for the covariance (23)) 

(Tr(exp(Ty^))) <Ci(27rTcr2)-i(sinh(wT))"i [ dxdycxp{-l-V2{{x,y)) 



'' ds 



-E 



G{{x,y),{x,y)) + G{(\s,cis) 



(68) 



<C2(27rTCT2)-i(sinh(wT))-i / dxdyeyiY>{-~TV2{x,y)){x^+y'^)~^ 



Eq.(68) shows that the growing random electromagnetic field improves localiza- 
tion. As an example we consider 



V2 = \xny\ 



(69) 



(the case a = 2 has been discussed by Simon [42]). The classical criterion for a 
discrete spectrum ( eq.(68) with 7 = 0) is not satisfied ( the region in the phase 
space with the classical energy less than E has an infinite volume, see [14], [42]). 
However, any 7 > (random vector field with a growing covariance) leads to a 
finite trace . Note that the results of [20]- [21] concerning the discrete spectrum 
do not apply directly to the vector potential (23)and the scalar potential (69) 
because the covariance of the magnetic field B is decaying as |x|^'''^^ (7 < 1). 
Hence, it is bounded in the mean. 

In the model (59) (with V3 = 0) we can obtain some estimates on the off- 
diagonal of the heat kernel as well. Let Kr{p\ x, y, ; x' , y') be the Fourier trans- 
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form of Kr{x, y, z; x', y' , z') in z' — z then 



Kr{p; X, y, ; x' , y') = (27raV)-i cxp(--i^(x - x'f - -^{y - y'f) 



2Ta2 



2Ta2 



E 



exp 



(^(p + A3(q(.)))2 + y2(q(.s)))d. 



< (27rCTV)"^ exp(-- — -(x - x'f - 



(70) 



2tct2 



2Ta2 



(y-y'f) 



^ds 



-£; 



exp 



2hm 



(p + As{q{s))f--V2{q{s)) 



The expectation value over Ag, can be calculated exactly. Let us consider a 
simple case of a translation invariant Gaussian field with G(x, x') = G'(x — 
x') (there is no scale invariance if G(0) is finite). After a calculation of the 
expectation value on the rhs of eq.(70) we obtain 



1 



{Kr{p; X, y, ; x', y')) < {2TTa^T) exp(-— ^(x - x'f - 



1 



2TCr2' 



2TCr- 



■iy-y'r) 



T 



exp ( - ip2(!^ + G(0))-i - J^2(q(.)))] (1 + ^G(O))-^ 



(71) 

When r ^ oo then the upper bound of eq.(71) is decreasing as exp(~ 2G('o'| P^) 
for a large p. We could interpret such a decay of the Fourier transform of the 
heat kernel as a confirmation of the behaviour exp(— ai?i~ (z — z'f) (which has 
also a Gaussian Fourier transform) in the lower bound (25) of Theorem 2 (scale 
invariance of the vector potential has not been assumed there). 

The decay of Green functions is important for correlation functions of the 
complex scalar fields interacting with an electromagnetic field in the Ginzburg- 
Landau model 

(0*(x)^(x'))-(a(^")(x,x')) (72) 

where Qm ' is defined in eq.(37). 

The lower and upper bounds on the higher order correlations of the scalar 
fields can be studied by means of our methods as well. In such a case the integral 
/ Adq must be extended to many paths joining the points x^ as the arguments 
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of the scalar fields </). For example 

(r(x)r(y)0(y')0(x')) = (g(^'")(x,x')a(^'°ny,y')) + (x - 



y) 



= / circir'(2^Ta2)-*(27rr'a2)-f exp(--l^(y - y')^) - :r^(x - x')^) 

{E exp (^- / A(qyyO o dq^yy + - A(qxx') o ^Qxx' j ) + (x ^ y) 

(73) 
where (x -^ y) means the same expression but with exchanged arguments. We 
can calculate the expectation value over the electromagnetic field and derive 
upper and lower bounds for the correlation functions (73). The important ques- 
tion to be answered is whether the decay of correlations holds true for any two 
points tending to infinity in the multi-point correlation functions of the scalar 
fields. This problem needs further investigation. 
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